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Propagation of waves in a compressible plasma, bounded by a rigid convex surface, is considered 
in this paper. The situation is idealized to the extent that Maxwell's equations, when combined 
with continuum theory of fluid dynamics, are separable. Specifically, the model is a perfectly conduct- 
ing cylinder of infinite length which is excited by a uniform voltage applied to an axial slot. It is 
shown that both electromagnetic and electroacoustic waves are excited in the plasma. Of particular 
interest is an azimuthal surface wave which circulates around the cylinder with exceptionally low 
attenuation. In the limiting case of a cold (incompressible) plasma, the surface wave is not excited. 

1. Introduction 

It is now known that certain significant differences exist between waves in cold plasma and 
waves in warm plasma [Stix, 1962; Denisse and Delcroix, 1963]. In the latter case, it is found, 
that in addition to the electromagnetic-type waves, one may also excite electroacoustic-type 
waves by a source. In homogeneous media, the propagation characteristics of these two wave 
types are determined bv uncoupled wave equations. The presence of either inhomogeneities 
or a d-c magnetic field will produce some coupling between the electromagnetic and electroacoustic 
waves. A particularly simple example of coupling occurs when a semi-infinite compressible 
plasma is bounded by a rigid dielectric with a plane interface [Hessel et ah, 1962; Seshadri, 1964; 
Wait, 1964a]. A plane electromagnetic wave incident from either side of the interface will produce 
secondary waves of both electroacoustic and electromagnetic types in the plasma. A notable 
exception occurs when the electric field of the incident wave is purely parallel to the plane interface, 
in which case no coupling occurs. 

An important parameter emerging in the analyses for electroacoustic-electromagnetic effects 
is the ratio (u/c) where u is the velocity of sound in the electron gas and c is the velocity of light 
in vacuum. Because this ratio is very small (e.g., 10 -3 or 10" 4 ), in both laboratory and ionospheric 
plasmas, the electroacoustic effects are usually negligible insofar as they modify the electromagnetic 
characteristics of the waves. However, there appear to be certain configurations where the electro- 
acoustic effects are particularly noticeable. Such an example considered in this paper is a 
compressible plasma medium bounded by a convex rigid surface. A source located on the surface 
will excite the usual creeping-type electromagnetic waves but, because of the compressibility, 
a strongly trapped surface wave will also be excited. On physical grounds, one may expect that 
deep in the "shadow" the trapped wave component would dominate the creeping wave components. 
This conjecture is verified by the analysis in the present paper. 

2. Formulation 

To simplify matters, a rather idealized model is considered. The d-c magnetic field is ne- 
glected. The convex surface is taken to be an infinite circular cylinder of radius a immersed in 
the plasma. The source is a narrow axial slot which is excited throughout the length of the cylinder 
by a voltage Vq. This configuration is almost completely analogous to the problem of a conducting 
sphere excited by an annular slot. Since the formal results for the sphere problem were given in 
some detail in an earlier paper [Wait, 1964b], the derivation of the corresponding formal equations 
for the cylinder problem will be omitted. The emphasis here will be on the transformation of the 
solution to a meaningful and possibly useful form. 
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The plasma medium is regarded as a one-component electron fluid. In other words, the ions 
are neglected in the equation of motion, yet their presence is required to neutralize the plasma. 
It is also assumed that the amplitude of the electromagnetic and the acoustic waves is sufficiently 
small that a linearized theory is valid [Oster, I960]. The average density of the particles is n and 
this is regarded as a constant in the plasma region. The pressure deviation of the electrons from 
the mean is p and their mean velocity is v. As usual, the electric and magnetic fields are denoted 

E and H, respectively. Collisions between particles and other forms of damping are consistently 
neglected in this paper. 

The linearized hydrodynamic equation of motion is 

dv ~£ <-, 

mno—=n eE—Vp, (1) 

where e and m are the charge and the mass of the electron, respectively. The equation of contiiv 
uity, when combined with the equation of state, leads readily to 

u 2 mn \/ - v=—dp/dt, (2) 

where, as mentioned, u is the velocity of sound in the electron gas. 

Maxwell's equations for the electromagnetic fields in the plasma are given by 

S7xE=-HodH[dt 9 (3) 

and 

V X // = eodE/dt + n ev, (4) 

where /x and e are the magnetic permeability and dielectric constant of free space, respectively. 
[Note that for electrons, e= — |e|] 



3. Formal Solution 

Without subsequent loss of generality, all field quantities are assumed to vary as exp (iwt) 
and, thus, the derivative d/dt may be replaced everywhere by ico. 

The conducting cylinder is now alined to be coaxial with the z-axis of a cylindrical coordinate 
system. Thus, the region p > a is the compressible plasma. Also, because of the two-dimensional 
nature of the problem, the derivative d/dz may be set equal to zero. 

With the simplifications introduced above, it is now a simple matter to show that the electric 
and velocity field components may be found from the scalar pressure p and the magnetic field 
which has only a z component. Thus, for p> a, 

F — - dffiz g dp 

p iewp d(f> eo) 2 m dp' '**' 

F = 1 8H Z e dp 

160) dp eco^mpdcj) W 

_ e dH z €p dp ^ 

p —€(jt) 2 mp d(f) ieconomdp ^ ' 

e dH z €o dp 

v <t> = — T""E ~ 7X' (8) 

eorra dp lecojiomp d<p v ; 

where — = 1 ; and cog = 
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The p and H z satisfy 

where ^ = /x 6w 2 and k 2 = {a)/u) 2 (e/eo). 

Solutions of (9) may be written in the form 

p = 2AnHM(k p p)e-™ ( io) 

n=-oo 

and 

ff*= J CJfi*Kkep)er«+ 9 (11) 

n = - x 

where H ( %\Z) is the Hankel function of the second kind of order /i with argument Z, while ^ n 
and C n are constants yet to be determined. The use of the Hankel function of the second kind 
assures outgoing waves at infinity. As indicated, the summations are over integer values of n 
so that periodicity in (/> is assured. 

On the assumption that the cylinder is acting as a rigid body, the boundary condition on the 
velocity is v p = at p = a and < </) < 2tt. Thus, using (7), (10), and (11), it follows that 

An_ n en H (2 ^{k e a) 



C n w€o (k p a)H^'(k p af (12) 

where the prime indicates a derivative with respect to the argument of the Hankel function. 
The coefficient A n may be found from the prescribed condition of the tangential electric field 
E<t>(a) on the cylinder. Explicitly, for p = a, 

E*=E<>(a) =^ for - (A/2) < <f> < (A/2) ( 1 3) 

= 0for (A/2) < |0| ^77, 

where A is the angular opening of the slot. 

A simple exercise in Fourier analysis gives the alternate form, for p = a, 

ti,Ka) 2ira£„ (nA/2) (14) 

for the whole interval — tt < (j) < tt. As indicated, this representation is appropriate for a uni- 
form electric field of strength Vo/(aA) within the slot. A more realistic variation [Wait, 1959] 
of the electric field within the slot is given by 

E * M = air[(A/ 2 ft-0»]* f ° r " A/2 < * < A/2 ' (15) 

where Vq is still the voltage across the slot. The corresponding Fourier representation for this 
variation is readily found to be 

£*(a)=-|2- J / (nA/2)e-** (16) 

n= — oo 
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where Jo is the Bessel function of the first kind of order zero. In general, we may write 

£*(a)=^ § /(nA/2)e-«, (17) 

where f{x) = (sin x)\x for the uniform field assumption and f{x) =Jo(x) for the nonuniform field 
in the slot. In both cases, f(x) approaches unity as nA/2 tends to zero. In this limiting case, 

£ * (a) = £2 e "^ = ^ 8W>) ' (18) 

n= — x 

where 8{4>) is the unit impulse function at = 0. In what follows, the function /(Aw/2) is retained 
in the analysis. 

Using (6), along with (10), (11), and (12), it is now a straightforward matter to apply the pre- 
scribed field condition at p = a. This process yields 

K»eea> /(nA/2) 

27Tk e aH\?Xkea)(l-8„)' (19) 



whe 



Thus, 



_ rcVoH 2 H»Kkpa) m\k e a) 
" (lea) (k p a) Hf'(kpa) Hf'(k e a) ' (20) 



Voieoj - /(nA/2) H<*\k e p) e -M 
"' 2TrkeaZ(l-8„)HW(kea) 



(21) 



and 



Kotcnoe ft »/(nA/2) H^(k„a) H^(k p p) 

P 2ne (k e a)(k p a) n ^ x (1 — 8„) Hf'(k e a)Hf'{k l> a) ' {ll) 

If the plasma were to become incompressible (i.e., cold), the parameter (k p a) would approach 
infinity. For this limiting case, the pressure p— >0 and the parameter 8 n ~ > 0. The resulting 
expression for H z is then identical to the magnetic field of an axially slotted cylinder in a dielectric 
medium of dielectric constant e and wave number k e [Wait, 1959] . 

Equations (21) and (22) are exact expressions for the magnetic field and the pressure, respec- 
tively, in the compressible plasma. Corresponding expressions for the electric and velocity field 
components may be obtained from (5) to (8). Some simplification is achieved in the far field when 
p — > oc. Then, the Hankel functions of arguments k e p and k p p may be approximated in the manner 

// " 2,(A(p) ~ (dip) ' 2f " exp (_ ikeP) < 23 ) 



and 



/ i \ 1/2 
H n ] (k P p) = \^f-J i n exp (- ik p p) . (24) 
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4. Application of Watson Transformation 

The subsequent discussion will deal with the behavior of the function H z . In particular, 
we shall develop an expansion for H z which should give some physical insight into the nature of 
the diffraction phenomena associated with the rigid cylinder. 

For convenience, we write 

H z = J /i„e-'"*, (25) 



hn 2nx HW(x) 1-S„ ' (26) 

with x — k r a and %\ — k e p. 

Following the prescription of the Watson transformation [Sommerfeld, 1949; Wait, 1959], 
the summation in (25) is replaced by a contour integral in the complex v plane, such that 

H z = — I h v e~^ -?— dv, (27) 

li J sin vtt ' 

G + C 2 

where the poles of the integrand are at v = n. Here, C\ is a straight line, just below the real axis, 
which runs from left to right; C2 is a straight line, just above the real axis, which runs from right 
to left. It is a simple matter to verify that 2m times the residues of the poles in (27) leads to the 
summation (25). Thus, the equivalence is established, provided that C1 + C2 enclosed only the 
poles at v — n. 

It is demonstrated below that h v has a number of complex poles which lie in the fourth quadrant 
at v = Vs where s indicates the number of the pole. Because of symmetry, poles also exist at 
v = —v s in the second quadrant. Furthermore, C\ may be closed by an infinite semicircle in the 
lower half plane of v without changing the value of the integral. Similarly, C> 2 may be closed by 
an infinite semicircle in the upper half plane. Thus, H z may now be replaced by a summation 
of the residues of the poles v 8 and — v g . Therefore, 



H z = -2n^ 



1 cos v h {tt — (\>) 



A-Jl 

dp h v 



(28) 



where the summation is over the poles v % which are solutions of 

(l/^) = 0. (29) 

The latter pole-determining equation is exactly equivalent to the equation 

8*=!. (30) 

where 

xy H^{x)H^\y)- (31) 
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where x = k e a and y = k p a. It is immediately evident that an evaluation of (31) leads to the propa- 
gation constant for the waves which "circulate" around the cylinder. Individually, they contain 
the factor exp (— iv s (j>) for waves traveling in the positive (f> direction or exp (-hi*/ g <£) for waves 
traveling in the negative </> direction. If is measured in a positive sense it is convenient to write 

cos v s {rr — (b) 

± St = fc-M&fo), (32) 

sin V S TT v ' 



where 



C s (4>) = 



l_|_ e -i2y s (7r-4>) 

1 + e' 21 ^ 



(33) 



It is evident that when <f) <tt and provided — Im v s tt > 1, the function G s may be replaced by unity. 
Under this condition, waves circulating around the cylinder are highly attenuated. 

A study of (30) indicates that there are two sets of poles. The first set which we shall describe 
as the electromagnetic or EM type, while the second set is called the acoustic or A type. The 
EM-type modes are located where \v\ is of the order of x, while the A- type modes are located where 
\v\ is of the order of y. Recognizing that y > x, it is evident that a different type of approximation 
is to be used in the two cases. With this point in mind it is desirable to rewrite (28) in the form 



H,=m+m, 



(34) 



where 



m 



-2ni2 



gfW 



■ - d 2-1 ex P ( — «1w. 

dv h v \ e 



(35a) 



and 



#«=-2tt;2 



GM>) 



I"- 1 

dv h v a 



exp (- ii/«0). 



(35b) 



Here, the superscript e refers to the contribution from the EM-type poles at v e s , whereas the super- 
script a refers to the contribution from the A-type poles at v*. 



5. Pole Determinations — EM Type 

For the EM-type poles, we may use the third-order approximation for the Hankel functions 
such that 



and 



i /2\ 1/3 
i /2\ 2/3 



(36) 



(37) 



where t — (v — x)(2lx) 1/3 , w(t) is the Airy function defined by 



w(t) = exp (- 2m/3X- Tri/3) 1 ^) [(2/3)(- 1) 3/2 ] , 
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and w'{t) = dw(t)/dt. This third-order approximation is valid when v is in the vicinity of x which 
itself is to be much larger than one [Wait, 1962] . Within this same domain the Hankel function 
H i2 \y) may be replaced by the first term of its asymptotic series. Thus, H {2) '{y) = — iH (2 \y). 
Using the above approximations for the various Hankel functions, it follows without difficulty that 

K= ^^\2) w'{t) 1-8,' (38) 

where 

. . . /o>oV v 2 (xVl* w{t) 

and 

Y = (2lxyl*k e Zo, z = p-a. 

Again, in order to simplify matters, the height z of the observer has been assumed to be small 
compared with the radius of curvature a. 

Using the third-order approximations described above, the residue series representation for 
EM-type waves is found to be 

Ht-H^^WtPP (40) 

where 

and t 8 are roots of the equation 

Co 2 v 2 /xVl 3 

w ' {t)+i ^\V ^ )=0 ' (42) 

and 

»°-m (§)"*-• 

The quantity H° z may be identified as the field of a "delta" slot at distance a<\) on an equivalent 
flat ground plane. 

Equation (40), which is a fairly general result, may be simplified somewhat when </> is not near 
77 and provided the width aA of the slot satisfies k e aA <1. In this case, both G^((/>) and /(z-f A/2) 
may be replaced by unity. Furthermore, under the assumption y > x, the dimensionless factor 
q e has a magnitude which is small compared with unity. Thus, the factor t s — q 2 — q e (2lx) 1/3 in 
(40) may be replaced by t s to within a good approximation. Thus, the simplified version of (40) 
reads 

„- ecoVp /2\ 1/3 w{t s -Y) . l*>\UUJL-\ 

H * = ~ \x) T7M^~ eXP [-«*-*(*/2) 1/8 ^]- (44) 

To determine the solution of (42), an iterative method may be adopted. For the first-order 
iteration, v 2 is replaced by x 2 in the factor multiplying w(t). Thus, (42) may be written 

w'(t)-q e w(t)=0, (45a) 
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where q e is given by (41). A solution of (45a) appropriate for small values of q e , is written 

t s = r s -\ — q e ~ 0/ x » q 2 e + terms containing high powers of q e , 

T s ATs) 

where r s are solutions of 

m;'(t)=0. (45b) 

Values of r g are well known from the theory of Airy functions [e.g., Fock, 1945]. Thus, 

n = 1.0188 exp (-J77/3), 

72 = 3.2482 exp (-jtt/3), 

73 = 4.8201 exp (-itt/3), 

etc. 

A higher order solution of (42) may be obtained by replacing v in (42) by x + (x/2) tf3 t s where t s 
is given by (45b). Then, a higher order solution 7* is obtained from 

w'(T)-QMT)=0, (46) 

where 



^feiWW/! 






1/3 



where, as usual, x = k e a and y—k p a. 
Thus, the higher order approximation is 

r.-*+£ft-2(J3i<B+. • • • 

If one wished, further iterations could be effected. However, as a practical matter, the zero-order 
solution given by (45b) should be quite adequate. 

In the limiting case of an incompressible (or cold) plasma where y = k p a^> °°, the parameter 
q e vanishes completely. The root t s occurring in (40) is then identical to r*. The residues series 
representation then has the same form as the groundwave field for radio propagation over a spheri- 
cal earth in an atmosphere of dielectric constant e[e.g., ch. V in Wait, 1962]. 

The preceding development indicates that the EM-type residue waves, circulating around 
the rigid cylinder, are only slightly modified by the compressibility of the ambient plasma. 



6. Pole Determination — A Type 

We now turn our attention to the contribution from the electro-acoustic or A-type poles. As 
indicated, the poles occur when v is of the order of y and since y > x, the Debye approximation 
[Sommerfeld, 1949] for the Hankel function fftf (x) is appropriate. Therefore, it is found that 



where x = k e a and x\ — k e p = k e (a -+■ Zo). 



(47) 
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In view of the large magnitude of v, (47) may be further approximated to 



Hf\x x ) _ x 



as exp 



■(x\—x) 



k e a 



exp 



Zo 
— v — 

a 



H^'(x) v 
Within the same approximation the function 8^, defined by (31), is represented by 



(48) 



6 "- U/ y^l 2) ' 



(y) 



(49) 



If it is assumed, provisionally, that the second order or Debye approximation is valid for H^\y), 
(49) may be simplified to 



6\, = + 



o>o\ 2 ^ 1_ 

a)) y[(Wy) 2 - 



l]l/2 



The pole condition 8 U = 1 then leads to the simple result that 



v = y[\ 



4\ -1/2 



a V or/ 



(50) 



This corresponds to an unattenuated surface wave propagation along the rigid surface with a 
propagation constant 

,>2\-l/2 



&p[l-(Ww) 4 ]- 1/2 or^(l+7§ 






This is precisely what one would expect for a compressible plasma medium bounded by a plane 
rigid surface. The applicability of the result to a curved surface rests on the validity of the Debye 
approximation for the Hankel function Hf\y). Certainly, when co<}/co 4< ^l, this situation would 
be violated. Therefore, it is necessary to use a third-order type approximation for H { *\y). 

Using the Airy function representations for the Hankel functions, as given by (36) and (37), 
it is now found that for the A-type poles, 



^ M 1/3 /cooV w(t) 



»-(D(f) 



CO/ (*)' (t) 



where 



(51) 



?=(y-y)(2/y)" 3 . 
The contribution from the A-type poles may now be written 

s l $ Q a 



where 



and 



y \2. 



m^\ ^_hea ( va zj> 

~v-v s 
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(52) 



(53) 



(54) 



Here, t s are solutions of 

w'(t)-q a w(i)=0 9 (55) 

where 

q a =(yl2yiH<ool<o)*. (56) 

The root of (55) of greatest interest is the one which has a surface wave character. Since q a is 
essentially real, it may be shown that [Wait, 1964c] the surface wave root of (55) is 

?o = Re? + iIm?o, (57) 

where 

lq a Kql 6Zq 7 a (58) 



and 



Im? ~ 2ql exp 



_4 3 _ 7 31__ 

Qa Uql 48rf ' * • J ' 



V a 12^ 48 9 « • • • J* (59) 

When q a becomes greater than about three, the imaginary part of ?o is negligible. Then, 

to = q 2 a , (60) 

which corresponds to 

i/g = y + (y/2) ^0 = y[l + ( Wco) 4 /2] . (61) 

Provided (co /w) 4 is sufficiently small, this result agrees with (50). The height-gain function H 
for this surface wave mode is given approximately by 

//o = exp [— vpoja] = exp (— k p z ) , (62) 

which indicates the strong trapping or "clinging" effect. 

When q a becomes small the attenuation, which is proportional to — Im ? , is appreciable. For 
example, if q a is 1.00, — Im t = 0.1551, and if q a is reduced to zero, — Im % approaches 0.5094. 
This corresponds to the cold plasma limit. 

There are other roots of (55) which are somewhat analogous to the conventional EM poles 
discussed in connection with (44). The height-gain functions for these modes are given by (53) 
with 5 = 1, 2, 3, ... . However, it would appear that for all practical purposes, contributions 
from these roots are negligible. 

In (52), the function / for the A-type poles is approximately represented by 

Av«AI2)=j{kpaM2) 

sin (&»aA/2) c . .- . .11 

— — Ti — a ins — l° r the uniformly excited slot 

(k p ak/2) J 

— Jo(k p aA/2) for the nonuniformly excited slot. 

Because k p a is a very large parameter, the function \f\ will be somewhat less than unity unless 
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aA, the slot width, was made exceptionally small. This points up the difficulty in exciting the 
electroacoustic-type waves by a finite source. 

In the preceding discussion, it has been assumed implicitly that a>o/cu 2 < 1. Thus, the propa- 
gation constants k e and kp are real. However, if a^/co 2 > 1, k e and k p become purely imaginary. 
The analysis given in this paper is applicable to this case if the factor (1 — w 2 /co 2 ) 1/2 is replaced 
everywhere by — &(a>J/(0 2 — 1) 1/2 . An interesting consequence is that the modes v\ and v a s have 
large imaginary parts except the surface wave mode i/g which is of the acoustic type. As seen 
by (50), this particular mode does not become highly attenuated as cojj/co 2 increases beyond unity. 

7. Concluding Remarks 

From the analysis given in this paper, it is evident that the finite compressibility of the plasma 
will modify the characteristics of the electromagnetic field in the vicinity of the rigid cylinder. In 
fact, for excitation by an axial slot of infinitesimal width, a strongly trapped surface wave is excited 
which travels around the cylinder with exceptionally low attenuation, provided the plasma is a 
non-dissipative medium. 1 However, the excitation of the surface waves is quite weak compared 
with the excitation of the electromagnetic-type modes. (As seen by (40) and (52), the ratio of the 
excitation factors is approximately (u/c) 213 where u is the velocity of sound waves in the electron 
fluid.) Nevertheless, if the field is measured within a distance A;" 1 from the surface, the field 
amplitude of the surface wave may be comparable with, or even much greater than that of the 
electromagnetic-type modes. Of course, in the limiting case of an incompressible (or cold) 
plasma, the surface wave is no longer excited. 



1 This means that collisionless damping is also neglected. 
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